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SHARP ANISOTROPIC ESTIMATES FOR THE BOLTZMANN 
COLLISION OPERATOR AND ITS ENTROPY PRODUCTION 

PHILIP T. GRESSMAN AND ROBERT M. STRAIN 



Abstract. This article provides sharp eonstructive upper and lower bound 
estimates for the Boltzinann collision operator with the full range of physical 
non cut-off collision kernels (7 > —n and s G (0,1)) in the trilinear L^(R") 
energy {Q{g, f), f)- These new estimates prove that, for a very general class 
of g{v), the global diffusive behavior (on /) in the energy space is that of the 

O^ 1 geometric fractional derivative semi-norm identified in the linearized context 

in our earlier works I15II16I . We further prove new global entropy production 

estimates with the same anisotropic semi- norm. This resolves the longstand- 

, ing, widespread heuristic conjecture about the sharp diffusive nature of the 

■^IJ ■ non cut-off Boltzmann collision operator in the energy space L^(R"). 
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1. Introduction and main results 
Our motivation for this study is derived from the physical Boltzmann equation, 

where the unknown f{t,x,v) is a nonnegativc function. For each time t > 0, 
the solution /(t, •, •) represents the empirical measure of particles. The spatial 
coordinates are generally x G fl and the velocities are v G M", where fl C K" 
is a domain (n > 2). The Boltzmann equation, derived in 1872, is one of the 
fundamental equations of mathematical physics and, in particular, a cornerstone of 
statistical physics. 
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2 P. T. GRESSMAN AND R. M. STRAIN 

In this article, we prove new sharp, anisotropic fractional derivative estimates 
for the Boltzmann collision operator Q. This is a bilinear operator which is given 

by 

(1) Qig^Div)"^ I dvj daB{v~v,,a)[g'J'-g*f]. 

We arc using the standard shorthand / ~ f{v), 5* = ff(v*), /' — f{v'), g'^ ~ 9{v'^)- 
In this expression, u, w* and v' , v'^^ are the velocities in R" of a pair of particles 
before and after collision. They are connected through the formulas 

(2) .' = !i±^ + t^., < = ^-^^a, a.S"-. 

These formulas correspond physically to elastic collisions with conserved quantities 



Our main focus is to study the sharp, anisotropic fractional diffusive effects induced 
by this operator under fully general physical assumptions on the collision kernel. 

The Boltzmann collision kernel _B(u — u», a) for a monatomic gas is a nonncgative 
function which only depends on the relative velocity \v — w*! and on the deviation 
angle 9 through cos0 = (k,a) where fc = (w — w,)/|w — w»| and (•, •) is the usual 
scalar product in R". Without loss of generality we may assume that B{v — v^,a) is 
supported on {k, a) > 0, i.e. < < ^. Otherwise we can reduce to this situation 
with the following customary "symmetrization" : 

B{v-v^,(t) = [B{v-v^,(t) +B('u-UH.,-o-)]l(fc^^)>o. 

Above and generally, 1a is the usual indicator function of the set A. 

The Collision Kernel. Our assumptions are as follows: 

• We suppose that B{v — w,, cr) takes product form in its arguments as 

B{v ~ u,,cr) = $(|u- v,|)6(cos6'). 

In general, both b and <!> are non-negative functions. 

• The angular function 1 1— >■ b{t) is not locally integrable; for C(, > it satisfies 

(3) ^<«in"-'^6(cos0)<-^, ,se(0,l), V0e(o,| 

Some authors use the notation v = 2s € (0,2), which is equivalent. 

• The kinetic factor z 1— >■ ^dzj) satisfies for some (7$ > 

(4) $(|iJ -u*|) = C$|w-iJ*p, 7 > -n. 

Our main physical motivation is derived from particles interacting according to 
a spherical intermolecular repulsive potential of the form 

0(r)=r-(P-i\ pe(2,+oo). 

For these potentials. Maxwell in 1866 showed that the kernel B can be estimated. 
In dimension n ^ 3, B satisfies the conditions above with J ~ {p — 5)/(p — 1) and 
s ~ l/(p— 1); see for instance p,[^. Thus the conditions in (jH]) and (|1]) include all 
of the potentials p > 2 in the physical dimension n = 3 as a particular case. Notice 
that the Boltzmann collision operator is not well defined for p = 2, see [25] . 

The celebrated Boltzmann _ff-theorem is one of the hallmarks of statistical 
physics. We define the H-functional by H{t) = — J^dx J^„ dv f log / for a suit- 
able domain Jl. Then the Boltzmann H-theorem predicts that, for solutions of the 
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Boltzmann equation, the entropy is increasing over time; formally (neglecting the 
boundary of fl), this corresponds to the statement 

This is a demonstration of the second law of thermodynamics. Here the entropy 
production functional, which is nonnegative, is given by 

(5) D{g,f)^^~[ dv Q{g J) log f. 

Moreover, this functional is zero if and only if it is operating on a Maxwellian equi- 
librium, e.g. (Hn]). This predicts that the Boltzmann equation exhibits irreversible 
dynamics and convergence to Maxwellian in large time. 

In our recent works [15l[T^ on the global-in-time stability of the Boltzmann 
equation with the physical collision kernels ^ and ^ , and near Maxwellian initial 
data, we introduced into the Boltzmann theory the following anisotropic norm: 

l/l^=.."'l/liv, + [ dv f dv' i{v) iv'))"^^ }/' 'Ss M^^^'Xi- 

T-+^= Jr" Je" d(w,w')^ 

Here we use one of the weighted L^ spaces, € G M, with norm given by 
mi^"^ I dv {vf Ifiv)]", l<p<oo. 



The weight is (v) = ^/T+\v^. If i? = we wiU write |/|lp = L/Ilp- We also record 
here the "dotted" semi-norm 



(6) l/U^,. = / dv dv {{v){v)) - ^ U[vy)<i- 



d{v,v'y 



The fractional differentiation effects are measured using the following anisotropic 
metric d{v, v') on the "lifted" paraboloid: 



d{v,v') = J\v - i;'|2 + i (|w|2 _ \v'\^y 



The inclusion of the quadratic difference \v\'^ — \v'\'^ is an essential component of 
the anisotropic fractional differentiation effects induced by the Boltzmann collision 
operator; it is not a lower order term as we will see in the following. Heuristically, 
this metric encodes the anisotropic changes in the power of the weight, which are 
non-locally entangled with the fractional differentiation effects. 

We have shown in [THl Section 2.3] that N'^^'' sharply characterizes the Dirichlet 
form of the linearized collision operator. In this work, we will show, perhaps more 
interestingly, that the anisotropic diffusive semi-norm (jS]) sharply characterizes the 
diffusive effects of the trilinear energy of the non-linear collision operator ([T]) under 
general conditions. We will furthermore prove that the diffusive effects of the 
entropy production functional ([5]) are also globally coercively controlled by (j6|) . To 
proceed further, we denote the L'^{E."') inner product by (•,•), which because of the 
context should not be confused with the scalar product in M". 
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It is well known that the Boltzniann collision operator ([T]) exhibits parametrized 
diffusion. For this reason we study the operator Qg , which is defined such that 



{Qi9J),h)= dv dv, daBiv-v,,a)[g'J'-gj]h 

= f dv [ diu [ da B{v -v,,a) gj {h' - h) = (/, Qgh) . 

This follows from the pre-post coUisional change of variables (w,w*,cr) — > (w',w^,fc), 
with unit Jacobian, which is standard |25j . We recall the widely used decomposition 
(from m, n) (/, Qgl) = -NgU) + Kgif) where 

NgU) - I [ dv [ dv, [ da B{v - v,,a) g, (/' - ff, 



Kgif) = ^ f dv f dvJ da B{v ~ v,,a) g, [(/')' - iff] ■ 

Now it is known from the cancellation lemma [T] that the second term Kg{f) does 
not differentiate at all, as will be seen in Section 11.31 below. 

In what follows, we make two assumptions on g = g{v) > from ([7]): 

Assumption U. For the upper bound inequalities below we suppose that g satisfies 
(8) f dv,\v~vX{v*y\gM\<Cg{v)% ae[j,^ + 2sl 3Cg > 0. 

JR" 

If s Cz (0, 1/2), we take i = 1, and if s G [1/2, 1), we use the power i = 2. 

Regarding Assumption U: if a > or if |u — w*| is replaced by the regularized 
kinetic factor {v — w,) then ([8]) will automatically hold with Cg « Ipj^i . Alter- 
natively, for any a > —n, condition (|S]) is satisfied, for example, by any bounded 
function which decays at infinity polynomially faster than order i + n (which will 
be the case whenever g belongs to any weighted Sobolev space of sufficiently high 
regularity and sufficiently rapid growth of the weight at infinity) . For the coercive 
lower bounds in Theorems [1] and [3] we assume: 

Assumption L. Let R > 6 > be fixed. Suppose the nonnegative, measurable 
function g satisfies 



(9) / dv, g* > Cg, 

■JBr\Ts 

uniformly for some positive constant Cg where Bfi is the Euclidean ball of radius 
R centered at the origin and Tg is any linear tube of radius d. 

Notice that such a 5 in Assumption L need not belong to L^(M") so long as it 
is locally integrable. In fact, one can show that any nonnegative, locally integrable 
g will satisfy ([9]) unless it equals zero almost everywhere on B^j. This is consider- 
ably more general than the assumptions used in [Tj , which prove previously-known 
entropy production estimates. It includes functions g which satisfy only "local 
conservation laws" for the Boltzmann equation as the following Corollary shows: 

Corollary L. Suppose g = g(v) > is a function belonging to the spaces L^{Bj^) 
and LlogL{Bii) with non-zero norms. Then ([9|) and (|lip hold for a constructive 
constant Cg > depending on R and the norms just mentioned. 
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Moreover, the Boltzmann H-theorem predicts that solutions to the Boltzmann 
equation converge as t — > oo to the foUowing Maxwelhan equihbrium states: 

(10) ^(p,^,T)(«)^=^'^^-A_exp^ 2^ 

where p, u, and T are the density, mean velocity and temperature of the gas re- 
spectively. Thus it can be expected that for p > and T < oo our assumptions 
(jH} and ([9]) will be satisfied for all time by any sufficiently regular solution of the 
Boltzmann equation which respects the H-thcorcm. 
We are now ready to state our main results. 



Theorem 1. (Main coercive inequality) Suppose that ^ holds. Then 



(11) N,{f) > C,\f\', 



where Ci ~ Cn,R,sCg /\g\Li(^Bji) CL^d Cn.R.s > only depends upon n, R, S, ^ and 
(|3|) . In particular, ^ and ^ together imply 

(12) -{Q{9J)J)>C^\f\%,,-C2C,\f\l., 

where all the positive constants, in particular Cn.B.s and C2, are constructive. 

We also have the following upper bound estimate: 

Theorem 2. (Main upper hound) Under (|8]) we have the upper hound estimate 

m{gJ),h)\<C:iC,\f\N^n\hU^... 

Here, as ahove, the constant C3 > is constructive as well. 

Collecting the estimates in Theorem [1] and Theorem [2] we obtain 

Ci|/|^.„ <-{f,QJ)+C2Cg\f\l. < iC, + Cs)Cg\f\%s,-,. 

These theorems thus show that the sharp global diffusive behavior of Qg in L^ is that 
of the geometric fractional semi- norm ^ , up to terms which do not differentiate. 
Furthermore we can coercively control the entropy production functional ([5|). 

Theorem 3. (Entropy production) Given ([8|) and ([9]), we have the constructive 
lower hound 

^(5,/)>Ci|//|^,,,-C2C3|/Ui. 
The constants Ci, C2 arc same as those appearing in (J12p . 

Each of the above anisotropic fractional derivative estimates improves upon 
previously-known estimates, which were formulated in various local and global 
isotropic Sobolev spaces such as those found in [T|[2|l4l[5l[T2 | [T9 l [24] . More discussion 
may be found in Section ll.il 

In what follows, we use the notation ^4 < B to mean that there exists a finite, 
positive constant C such that A < CB holds uniformly over the functions and sum- 
mation indices which are present in the inequality (and that the precise magnitude 
of the constant is unimportant). The notation i? > A is equivalent to A< B, and 
A^i B means that both A<B and B < A. 
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1.1. Historical remarks. Cercignani [H p. 85] in 1969 (about forty years ago) 
noticed that the linearized Boltzmann collision operator L in the case of the Maxwell 
molecules collision kernel {p = 5) behaves like a fractional diffusive operator; more 
precisely, in [5] it was noticed that the eigenvalues of L grow like —rn^i^ as tti — >■ oo. 
Over time, this point of view transformed into the following widespread heuristic 
conjecture on the diffusive behavior of the Boltzmann collision operator ([T|): 

./ >-> Q(ff, /) ^ —{—AyYf + lower order terms. 

Note that (—At,)* is a flat fractional Laplacian. See for example [l]|5l[25j. This 
point of view is well-known to be correct locally [1] in, for instance, a ball Bfj C K" 
with < i? < oo. Our main motivation is to provide a global, sharp anisotropic 
correction to this heuristic conjecture in terms of energy estimates for the trilinear 
energy under ([8]) and ([9]). We furthermore believe that the information provided 
by these new estimates in Theorems [l] [21 and [3] will be quite useful to future work 
in a wide variety of contexts within the Boltzmann theory. 

We also point out that the best possible comparisons of the space | • In"--! to the 
weighed isotropic Sobolev spaces Hg{R^) are established by the inequalities 

(13) I'^li^+.^CK") + l^lff^(K") ^ l^^l^- ^ l^lff; + ..(K")- 

Here i//(R") '=' {/ e LJiR-) : |/|2^.(«„) ^=*' 4„ dv {vf \{I - A..r'V{v)\^ < ^} 
is the standard isotropic fractional Sobolev space. These inequalities can be easily 
established as in [161 Eq. (2.15)]. 

Historically, anisotropic behavior in the linearized context was noticed as early 
as Pao [21] in 1974, where he studied the symbol of the Fourier transform of the 
linearized collision operator. Indeed, a key difficulty in the analysis of [H] was the 
presence of the cross product of the frequency (derivative) variable and the velocity 
(weight) variable; Pao's pseduodifferential operator approach required an intricate 
knowledge of Bessel functions. Further studies of the Fourier transform of the Boltz- 
mann collision operator were given, for example, in [I1I31I71I12] and the references 
therein. The sharp anisotropic differentiation effects for the Dirichlet form of the 
linearized Landau collision operator have been studied, for example, by Guo [17| in 
2002, and Mouhot-Strain [20| in 2007. The Landau operator involves whole deriva- 
tives rather than non-local, geometric fractional derivatives. It is worth discussing 
briefly the terms "anisotropic" and "isotropic" because, for example, Pao [21] uses 
the term "isotropic" to describe the effects that we are now calling "anisotropic." 
Additionally, the terms anisotropic or non-isotropic are sometimes used loosely to 
mean there is a gain of weights or a gain of regularity in the relevant estimates, 
which does also occur in some isotropic situations. We use the term anisotropic in 
this paper to mean that there is an essential coupling between the non-local "direc- 
tion of differentiation" and the power of the velocity weight at infinity (meaning, for 
example, that the spaces H^ and ^^+2s appearing in p^ are sharp but not equal 
to one another) . This strict notion of anisotropy can be observed in the Boltzmann 
theory using delicate calculations involving the Fourier transform; see, for example, 
[511^ . It is worth noting that a large number of alternate formulas can be derived 
for the Fourier transform of the Boltzmann collision operator. However since the 
Fourier transform is a Euclidean invariant object, it seems to be difhcult to use the 
Fourier transform point-of-view to prove sharp anisotropic energy estimates in the 
presence of essentially non-Euclidean geometries. 
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In this paper, we treat the Bohzmann colhsioii operator as a fractional, geometric 
Laplacian with the geometry of a "Ufted" paraboloid in R"+^ which was introduced 
in our previous work [15[|16| . In particular, we do not use the Fourier transform 
at all. With this point of view, we can apply a generalization of Littlewood-Paley 
theory, as in Section |3l to prove our sharp anisotropic upper bound estimates in 
Theorem[2l Our Littlewood-Paley projections, developed in [161 . arc inspired by the 
Littlewood-Paley-Stein theory of Stein [22], the geometric Littlewood-Paley theory 
of Klainerman-Rodnianski [TH], and the original physics representations of (IT|) in 
terms of delta functions on the collisional conservation laws. 

Recently this point of view has played a crucial role in our proof [T5l[T6l[23] of 
global existence of unique, non-negative, near-Maxwellian classical solutions to the 
Boltzmann equation (1872) which exhibit rapid convergence to equilibrium. These 
results cover the full range of physical collision kernels derived by Maxwell in 1866 
from an inverse power intermolecular potential, and they additionally resolve a 
conjecture from [20]. Notice also [6j for other related results. 

The sharpest previously-known estimates for the collision operator arc expressed 
in terms of isotropic Sobolev spaces and correspond to the inequalities 

-{Q{9J)J) > c(g)|/|^ c|5b ^^^ ^,m\f\l\im- 

Here (7-|-2s)"'" = max{0,7-|-2s}. For estimates of this type, see, for example, [TJI^J 
[5l ll2[[24l[25] and the references therein; in particular, these (and other) estimates in 
terms of isotropic Sobolev spaces are stated in [5] Theorem 2.1, Theorem 2.6]. Here 
C{g) > depends upon the L^r^ 2—y+}(^^) ^^^ LlogL(R'^) norms oi g and also 
the collision kernel B. These bounds are proved in [5j when <I> in (|4]) is replaced 
by the regularized kinetic factor <&(|2:|) = C^ {z^ which removes the analytical 
difficulties at zero and infinity. Notice that the weights appearing in these isotropic 
upper and lower bounds correspond to the best isotropic approximations of N'^''^ 
from above and below, as indicated by (|13p . Our estimates in Theorems [T] and 
[2] then improve upon these isotropic estimates in three ways. Firstly, we include 
the physical kinetic factors from ([3]) and ([4]). Secondly, we include the subtle but 
important sharp, anisotropic fractional diffusive effects in both the upper and the 
lower bound. And third, we only need the local information on g such as that in 
([9]) to prove the lower bounds. 

Regarding Theorem [S] many works study entropy production estimates in the 
non cut-off regime, as in, for example, [l][T2j[T9l[24] . These estimates have found 
widespread utility for both the spatially homogeneous and inhomogencous Boltz- 
mann equation, see e.g. [TI IT^I^ . The most widely- used estimate from [1] (which 
is sharp in comparison to ours locally) is the following: 

Dig, /) > Cg^R \^/f\%s(^BR) ~ ^°^^^ ^'^^^'^ terms. 

The lower order terms above can be found in [TJ Corollary 2]. Notice that the 
anisotropic effects of ([6]) are only present near infinity; our estimate in Theorem 
[3] thus implies this local estimate under more general conditions on g, as in ([5]), 
than those used in [T]. Theorem [3] is, moreover, a stronger anisotropic and global 
version of this local smoothing estimate [T] (stronger, that is, in terms of the global 
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anisotropic weight coupled to the fractional differentiation). We would like to re- 
mark that after the submission of the present results, we learned that some isotropic 
trilinear cocrcivity estimates were also obtained by Chen and He in |101lllj . 

Due to length constraints, it is difficult to provide an exhaustive set of references. 
However we refer to further references in our earlier articles [15[|16] and the reviews 
for more historical discussions of previous results. 



1.2. Outline of the article. The remainder of this article is organized as follows. 
In Section 11.31 we reformulate the Boltzmann collision operator ([T]) , using (O to 
reduce ([T2l) to ([TT]) : we also prove Corollary L. Then in Section[2]we begin to prove 
Theorem [2| we perform the anisotropic dyadic decomposition of the singularity and 
subsequently prove the "trivial" estimates as well as the cancellation estimates for 
the decomposed pieces. Following that, in Section |3] we describe the generalized 
anisotropic Littlewood-Paley projections which were developed in our earlier works 
[ISIITB] . We then perform the triple sum estimates for the trilinear form (p6| , which 
completes the proof of Theorem [2] Section |4] is devoted to the proof of (fTTj) . The 
methods used here are distinct from those appearing in our earlier works |15 1 ll6 j . 
In particular, we introduce new, non-oscillatory methods to derive a new inequality 
for Ng{f) which reduces the problem to the study of a family of convolution-like 
estimates. The key idea is to observe that Ng{f) and the semi-norm N'^''^ are 
already in a functionally similar form, but that the support of the integral appearing 
in Ng{f) is singular when compared to that of the semi- norm. By means of the 
convolution-like operation, we are able to "smear-out" the support of integration in 
Ng{f) from ([7]) and thereby make a direct, pointwise comparison. Then in Section[5] 
we prove Theorem lousing (|lip and a well-known splitting from, for instance, [71124) . 
Lastly, in the Appendix we derive the "dual formulation" (Q*f, h') for the trilinear 
form using a Carleman-type representation and methods from [TBI Appendix]. We 
will also derive a change of variables that we call the "co-plane identity" which is 
used in Section |4?T] 



1.3. The lower bound estimate p^ and Corollary L. From ([7]) and the 
cancellation lemma [1] Lemma 1] we quickly have the identity 

Kgif)^ f dv, g, {f * S){v,) ^ C I dv f I dv,g,\v-v,\^>Q, 



where (/^ * S) denotes convolution. In the formula above. 



(14) S{z) = ^4^ \zp r ' d9 sin"-2 e &(cos 9) 

2 Jo 



cos-(^+") - - 1 
2 



where now C" = ^^ J^^^ dO sin" "^ 61 6(cos6l) [cos-('^+") | - l] > 0. Note that 
C" is finite by virtue of ([3]), where Cn > is a dimensional constant. By ([8|) it 
follows immediately that Kg{f) < C" Cg 1/1^2 • Thus the term Ng{f) must contain 
all of the global, geometric fractional differentiation effects. Subject to ([TT|) . then, 
we have shown (fT2|) from Theorem [T] We will prove (fTTj) in Section H) 
Now we will give the proof of Corollary L: 
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Proof of Corollary L. Fix R > and let Tg be any tube of radius 6 > 0. By 
Jensen's inequality, we obtain 

meaning that 



\Ib dv^ g^\n{e + g^) I I \BnnTs\ J^ dv^ g^\n{e 

If iBfiCiTsl is sufficiently small relative to J^ dv^, g^ ln(e + .9*), then this inequality 
forces Jg ^rp 5* to be arbitrarily small with respect to /g du* g, ln(e + 5*). In 
particular, for small enough 6 (keeping R fixed), one can ensure that no more than 
half of the mass of 5* on the ball of radius R is contained in any tube Ts: 



/ dv^, gt> - du* g*, 

JBf,\T/i ^ Jb„ 



IBr\Ts ^ JBn 

for any tube Tg of sufficiently small radius. D 

2. Physical decomposition and individual estimates 

In what follows we prove all of our estimates for functions in the Schwartz space, 
5(M"), which is the well-known space of real- valued C°°(M") functions all of whose 
derivatives decay at infinity faster than the reciprocal of any polynomial. The 
Schwartz functions are dense in the anisotropic space N^''^; the proof of this fact is 
easily reduced to the analogous one for Euclidean Sobolev spaces by means of the 
partition of unity as constructed in [HI Section 7.3]. Moreover, in our estimates, the 
constants will not rely on the regularity of the functions that we are estimating. 
Thus, using standard density arguments, our estimates apply to any function in 
A^^'T or whatever the function space happens to be for a given estimate. 

We will now introduce the anisotropic dyadic decomposition of the singularity 
([3]) in Section 12.11 and all of the decomposed pieces of the relevant trilinear energy. 
Then in Sect ion [2T2] we perform the "trivial" estimates of the individual decomposed 
pieces, and in Section [2.31 we present the cancellation estimates. 



2.1. Dyadic decomposition of the singularity. Let {xk}kL~oo ^^ * partition 
of unity on (0, 00) such that |x/c|l~ < 1 and supp(x/£) C [2^*^^^, 2"*^]. For each k: 

Bu = B,{v - v,,a) "^^ mv -v4)b f/^--^a\] Xk{\v - v'\). 

Notice that we have the expansion 

(15) \v-v'\^ = J -— L 1- ( -,a) = |w-z;*psni^ 



\v-v,y /J ' ^' 2' 

2"^^ is equivalent to the condition that the angle 

_ ia r'/^mT-\OTaT-^l 



between a and "_^*| is comparable to 2 '^|w — u^j ^ . With this partition, we define 



(16) 



vXif,h)= f dv [ dv, [ da Bkiv-v,,a) gjh\ 
'D'l{f,h) = / dv dv^ / da Bk{v ~ v^,a) gt,fh. 
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Notice that both T)'^ and I?^ depend on g in spite of the fact that the notation 
used here suppresses this dependence. 

We will also express the collision operator ([T|) using the dual representation, 
denoted by Q* This is derived via a Carleman-type representation in the Appendix 
using methods from [ini Appendix]. With (pO]) in the Appendix we can see that 

(17) (/, Qgh) = (Q;/, h) "=^* / dv' f dv, f dn, B g, h' (/ - /') + 0,(/, h), 



Blv^V^ 2V-^,-^,. 



^dot„„_i^V "*'\2v'~v~v 



where the kernel B is given by 

(18) B 

\v' -v^\ \v- u,|"-^ 

and the operator O, above does not differentiate at all: 

(19) OM, h) '=' f dv' fh' f dv^ 5* / dn^ B f ^ ~ ^*j|"' ~ "f - 1 

Jm" Jr" Je-1 \ <P{v - v^)\v - v^\" 

We decompose this dual formulation as follows: 

Vlifji) = f dv' f dv, f dir, Bk g*h'f, 

JM" JR" J_E"' 

(20) /* 

'D'if.h) = / dv' / dv, / dn., Bk g.h'f, 

JR" JR" J Ei;'^ 

where we use the notation 

B (v - v.. ^^'-''-^'* 



■Bfe = 2 — -) iTTZ^Xfedi'-t' I) 

\v' — w*| \v — v,\"- ^ 

In the integrals (HZ]), ([T^ . and (PH)) . d-K^, is the Lcbcsgue measure on the (n — 1)- 
dimensional plane E^^ passing through v' with normal v' — v., , and v is the variable 
of integration (i.e., i?^ contains all v £ M" such that {v' — v, v' — «*) = 0). Notice 
the partition of unity guarantees that the kernel is locally integrable for any k. 

When f,g,h £ 5(M"), the pre-post coUisional change of variables and the dual 
representation p?|) . ((60)) yield the identities 



{f,Qsh) = {Q{gJ),h)^ J2 {V'iif,h)-Vlif,h)} 

k—~oo 

oo 

= a(/,/i)+ Yl {'Di{f,h)-v',{f,h)}. 

These will be the general quantities that we estimate in the following sections. 

2.2. "Trivial" estimates of the decomposed pieces. The next step is to esti- 
mate each of 2?^ , V'i , 2?^ and O, using only the known constraints on the size and 
support of Bk and Bk- This is what we will do now. 

Proposition 1. Under the assumption ([5|), for any integer k, we have the uniform 
estimate 

(21) \v'l{f,h)\<C,2'^'^\f\^.\h\^. 
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Proof. Given the size estimates for 6(cos 9) in ([3]) and the support of Xk , clearly 
(22) / d<j Bk < m\v - v,\) d0 0-i-2s <22s/c|^_^^p+2s^ 

Now we use the Cauchy-Schwartz inequality to obtain 



(23) \v'li.f,h)\ <2'^''' f dv [ dv, \v - v,['+'^' |<?*| \hf\ 

<2^^^([ dv]!]-" f dv, \g,\\v-v,\-r+A 



1/2 

dv |/ip / dv, \g,\ \v~v,\^+'^' 

Therefore ^ follows from dH). D 

Proposition 2. Given ^, the inequality below is uniform for any integer k: 
(24) \Vl{f,h)\<Cg2'-^'^\f\^.^^Jh\^.^^^^. 

Proof. As in the previous proposition, the crucial point in this inequality is the 
symmetry between h and / combined with Cauchy-Schwartz. This time the dual 
representation (|20p will be used. In this case, the key quantity is 

1 [ ^^^^(W-v.?-\v-vf\xk{\v-v'\) 



The argument of h in this expression follows from the identity 

/ V -V, 2v' -V - V, \ _ \v' - u*p - \v - u'p 
^ -' \|i;-w*|' |2i;'-t.-i;,|/ " \v - v'\^ + \v' - v,\^' 

The support condition yields \v — v'\ « 2"*"'. Moreover, since 6(cos6') vanishes for 
6 G [tt/2,tt], we have \v' — v,\ > \v' — v\. Consequently, the condition (|3]) gives 



\v' — v,\^ + \v ~ v'\^ J ^ Vb'~^*P, 

Notice that on E^^ we have \v — w*p = |u — v'\^ + \v' — u,p > \v' — w*p. Thus, the 
integral is bounded by a uniform constant times 

I,/ _ „ |n-l+2s 

d^^ I ;|Li+2. 1^' - y*r^'xki\v - v'\) < 2^^'\v' - v,\'^. 

By these estimates, it follows that 
(26) / d^„ Bfc<2'">'-^;*|T+2^ 

As a result, we obtain the upper bound 

\vtif,h)\<2'^'' f dv' f dv, \v'-v,\''+^^ \g,h'n 

JR" JR" 

Now the estimate ([M)) easily follows from (|5]) just as was accomplished in (P^ . D 
Proposition 3. Under the assumption ([5]) wii/i a — j, we have 

\0.{Lh)\<Cg\f\^\h\^. 
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Proof. We recall 0^,{f,h) from (fTQ]) . The key quantity here to estimate is the 
integral J^^, dny Bk{A-l) , where A "^^ ^^^("^I^;;!}^."" | L . Observe that, on £'^^ , 

\v'-v4' 



A = 
Now for any fixed a > 0, one has |c" — 1| < |c — 1| uniformly for < c < 1; thus 



/ dlTy Bfe|^-l| < / dlTy Bk 



\v-v'\'^ 
\v — u'p + \v' — w,p 



< / d^„ Bk min<jl, ' ' 



<min{l,2-2'=|z;'-t;,|-2} I d7r„ ^fe <2(2-'-»)fc|w'-w,|')'+2«-*, 

for both i = and i — 2 (where the last estimate in the series above follows from 
(El))- We conclude 






|7+2s 

fc: 2'=|i''-'(;^|<l 

fe: 2'=|li'-t)»|>l 

Now we complete the estimate with Cauchy- Schwartz and ([8]) as in (|23|) . D 



Proposition 4. For any integer k, under the assumption ^, we have the uniform 
estimate 

(27) |2?i(/»|<C, 2^^^ I/Il!., J/^IlL 



-.fc 



Proof. Notice that the operator 'D'^{f,h) is given by either ((20)) (with Carleman 
variables) or (fT6)) (without Carleman variables). By Cauchy- Schwartz, we have the 
inequality 

\'D'lif,h)\ < f dv f dv, f da Bkiv-v„a) \g,\ \fh'\ < ^f^M 
where Ai = /jj„ dw |/p /jj„ dw, /g„_i dcr -Bfc('y — v^,u) \g^\ and 
Ai = \ dv \h'\ / dw* / da Bk{v — v^,,a) jg*! 

dw |/ip / dw* / da Bk{v~v^,a) \g'^\ 

JR" JS"-1 

dw* |/i»| / dti / da Bk{v — v^^a) \g'\. 



The various equivalent formulas for A2 follow from the pre-post collisional change 
of variables (u, u,,cr) — > (u',u||,,fc), and the v <-> u», ct -H- ~a symmetry. Clearly 
Ai < Cg2'^'^\f\l:, as in (H]) and (l22|). For A2, we claim that 



7 + 2s 



(28) / d« / da Bk{v~v,,a)\g'\<2^''' dv' \g' \\v, ^ v' \ 



|7+2s 
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Then A2 < Cg2'^''''\h\'i2 follows from ^ and the proof will be complete. With 



7 + 2s 



the Carleman representation, [TBI Appendix] , the left-hand side of ([28|) is 
dv da Bk{v - v^,a) \g'\ = j dv' j diTy Bk \g'\. 

Now the claim ([25)1 follows from the estimate (pS)) . D 

2.3. Cancellations. In this subsection we study the differences (T>^ — P*^) and 
(P^ — V^Y The estimates herein will exploit the cancellations to obtain better 
dependence on k than in Scction [2.2l The price to be paid is that we must measure 
the magnitude of the differences anisotropically. 

The scaling dictated by the problem is that of the paraboloid: namely, that the 
function f{v) should be thought of as the restriction of some "lifted" function F of 
71 + 1 variables to the paraboloid (w, i|wp). Consequently, the correct metric to use 
in measuring the length of vectors in R" will be the metric on the paraboloid in 
]^n+i ijKiJuced by the (n+l)-dimensional Euclidean metric. To simplify calculations, 
we work directly with the function F rather than / and take its {n+ l)-dimensional 
derivatives in the usual Euclidean metric. 

To begin, wc will write down a formula relating differences of F at nearby points 
on the paraboloid to derivatives of i^ as a function of n + 1 variables. To this end, 
fix any two v, v' € R", and consider C : [0, 1] -^ R" and C : [0, 1] -^ R"+^ given by 

C(i9) = ^v' + (1 - d)v, and C(^) = i-dv' + (1 - d)v, ]- \^v' + (1 - ^)v\^\ . 

Here C lies in the paraboloid {(wi, . . . ,Wri+i) G R"^"^ | Vn+i = ^(t'l + • • • + t'^) }; 
also note that C(0) = v and C(l) — v' . Elementary calculations show that 



%^)^{v'^v,{a^),v'-v)), and ^ = {0,\v'-v 



2\ 



Now we use the standard trick of writing the difference of F at two different points 
in terms of an integral of a derivative (in this case the integral is along the path C) : 

(29) =l\^(^^.(^^F)im) 

where the dot product on the right-hand side is the usual Euclidean inner-product 
on ]R"+^ and V is the {n + l)-dimensional gradient of F. For convenience we define 



|VP-F(wi,. . . ,w„+i) = max sup 



^■vV F{vi,...,Vn+l) 



1,2, 



where ^ e R"+-'^ and |^| is the usual Euclidean length. 

If V and v' are related by the collision geometry, ([2]), then (w — u', u' — w,) =0, 
which yields that 

(C(i?), w' - w) = (w,, w' -v)-{l-d)\v- w'p. 

Thus, whenever |w — w'| < 1. which holds near the singularity (A: > 0), we have 



d'& 



< 
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In particular, for differences related by the collision geometry we have: 



(30) 



F v' 



- F V 



2 



<{v,)\v-v'\ / d^ |V|i^(CW). 

"'0 



Furthermore, by subtracting the linear term from both sides of ([29| and using the 
integration trick iteratively on the right-hand side of ((29|). we obtain 



(31) 



F v' 



,/|2 



F V 



dd 



(0) • VF{v) 



< {v.yiv-v'l''' I dd |V|''F(C(^)). 





We note that, by symmetry, the same result holds when the roles of v and v' are 
reversed (which only changes the curve C, by reversing the parametrization: (.{-d) 
becomes ^(1 — i^)). We will use these two basic cancellation inequalities to prove 
the cancellation estimates for the trilinear form in the following propositions. 

Proposition 5. Suppose h is a Schwartz function on K" given by the restriction 
of some Schwartz function H on R""*"^ to the paraboloid {v, h\v\'^). Let |V|*/i be the 
restriction of |VpiJ to the same paraboloid {i = 1, 2). Then, for any k>0, 



(32) 



\{V^+-V'l){f,h)\<Cg2 



(2s-i)k 



\I\Ll^JMhU 



Here when s € (0, 1/2) in ^ then i = 1 and when s € [1/2, 1) we have i = 2. 
Proof. For s G [1/2, 1), we write out the relevant difference into two terms 



h' -h 



h' 



d£ 
dd 



io)-ivh){v: 



di9 



iO)-iVh)iv) 



We further split {V^ - Vt) if,h) = V^ + V^^ where V^ corresponds to the first 
term in the splitting above. We begin by considering the last term 



pi 



dv 



du* 



dC 
da Bk{v - v,,a) gj -=(0) • {Vh){v). 
gri-i av 



dC , 



Notice that j5(0) is linear in v' — v and has no other dependence on v'. Thus 
the symmetry of Bk with respect to a around the direction "~^*| forces all com- 
ponents of v' — V to vanish (when integrated in a) except the component in the 
symmetry direction. Thus, one may replace v' —v with y^Z^' 



V, 



we obtain that 



= 0, the vector further reduces to " ^*i t^ — — 



^' I.-..I/- Since 
Since W~v\ ^ 2~^ 



V — v^, \v' — w|^ 



\v — f* 



dC, 



<2-"''>-w,|- 



The last coordinate direction of ^(0) is given by (w, v' — v) which reduces to 
'v, ---*\f~-\" \ <{\v'~v\' + \v- v.nv' - v\' (.,)) . 
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With these bounds for T)^^, we must control the integral 



pM <2- 



'2k 



dv I dih. 



da Bk 



V\h){l + \v-iu 



<2 



-2fe 



dv\ff / dv.,{v,)\g.A{l + \v-v.A-') 



da Bi 



S"- 



1/2 



dv{\V\h)' 



dv^{v^)\g^\ {l + \v-v^\ 1) 



da Bi 



S"-i 



1/2 



We complete the estimate for this term with ([22| and ([8]) for a = 7 + 2s and 
a = 7 + 2s - 1. Here we used that 2s - 1 > since s g [1/2, 1). (Note that ^ for 
the case a = 7 + 2s — 1 easily follows from the cases a = 7 + 2s and = 7.) 
We will now estimate I?". We use the difference estimate (|5T|) to obtain that 

|pi| <2-2fc f d^ f dv [ dv, [ da Bk {v,f \gj\ |V|2/i(C(t?)). 

Note that the factor 2^^'^ comes directly from ([3T|) . With that last estimate, 
Cauchy-Schwartz (as in the previous case), and ([8]), it suffices to show that 



(33) 



di9 



dv 



dv. 



da Bk{v,f\g,\ |V|2/i(C(^)) 



<Cg2'"=||Vr/i|j 



This bound follows from the change of variables u = 'dv' + (1 — 'd)v, which sends v 
to u. From the collisional variables ([2]), we see (with 5ij the Kronecker delta) that 



dui 

dVj 

with the unit vector k = 
du, 



{\-^)6^J+'^ 



dvj 



= 1- 



/c j ai , 



dvi 



{v - v,)/\v - v,\ 



-'^A 



2) '- + 2 
Thus the Jacobian is 



1 

2 



f (fc,^) 



Since b{{k,a)) = when (fc, cr) < from ^, and ■& E [0,1], it follows that the 
Jacobian is bounded from below on the support of the integral ([551 . But after 
this change of variables, the old pole k ~ (v — v,)/\v — v,\ moves with the angle 
a. However, when one takes k = (u — v,)/\u — w*|, then 1 — (fc, cr) « 1 — (fc, cr) , 
meaning that the angle to the pole is comparable to the angle to k (which does not 
vary with a). Thus the estimate analogous to (P^ will continue to hold after this 
change of variables, giving precisely the estimate in (j33|) . 

It remains to prove ([32)) for s e (0, 1/2). This estimate is exactly the same as 
the one for P" except that the cancellation term 7^(0) • {\/h){v) is unnecessary and 



we can use (|30p instead of ([?!]) which allows us to take i = 1 in ([5]). 



D 



Proposition 6. As in Proposition^^ suppose f is a Schwartz function on M" which 
is given by the restriction of some Schwartz function in M."''^^ to the paraboloid 
{v, ^l^l^) and define |V|*/ analogously for i = 1, 2. Then for any k > 0, we have 

(34) ICDl^V:) if,h)\<C, 2(2—)^- ||V|7Il^^,J/i|l^^,^. 

Here when s G (0, 1/2) in ([3]) then i = I and when s G [1/2, 1) we have i = 2. 
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Proof. This proof follows the pattern that is now well-established. The new feature 
in (p4|) is that, from ([20]) . the pointwise differences to examine are 

/-/'=(/-/'- §(1) • (V/)(«')) + ^(1) • (v/)K)- 

We again split (P'ji — I? J) (/, h) = V\+I)^J-, where V], corresponds to the first term 
in the splitting above. For the last term Pj,", we have 

pii ^=^f / dv' f dv, f dn^ B, g, h' ^(1) • (V/)(i/) = 0. 

To see this, note that, in this integral, as v varies on circles of constant distance to 
?;', the entire integrand is constant except for ^(1). If wc write 35 (1) as a sum of 
two vectors, one lying in the span of the first n directions and the second pointing in 
the last direction, it follows that we may replace the former vector by its projection 
onto the direction determined by v' — v^,. But since the original vector points in the 
direction v' — v, the projection vanishes. Since the last direction of ^(1) is exactly 
{v',v' — v), the corresponding integral of this over v also vanishes by symmetry. 
To estimate Vl we use ([30]) or ((3T|) to observe that 

|2?!| <2-*'= f d^ f dv' f dv, f dn, Bk {v,y\gM\ |V|7(CW) ■ 
Jo Jr^ Jr" J e^'. 

In other words, when s £ (0, 1/2), we use (PO]) and obtain i = 1 and, alternatively, 
when s £ [1/2, 1) in ([3|), we use (|3ip and obtain i = 2 above. Since V^^ = adding 
the cancellation term ^(1) • (V/)(w') when s £ (0,1/2) causes no new problems. 

Now we can estimate \'D\\ above using Cauchy-Schwartz, ([25)) . and ^ to conclude 

1/2 
that it is uniformly bounded above by a fixed constant times Cg multiplied by 

/ 1 \ ^/^ 

2^'-*^'>Il^, / dv,{v,y\g,\ j d§ f dv' [ diT, Bk |V|7(CW) ) • 
^+'= \Jr" Jo Jk" Je-1 J 

The integrals involving |Vp/(C('?)) are estimated as in (|33l) after changing from the 
Carleman representation to the a representation as in [I6j Appendix] . D 

3. Triple sum estimates for the trilinear form 

The point of this section is to prove Theorem [51 To proceed, we first recall 
the necessary anisotropic Littlewood-Paley theory adapted to the geometry of the 
paraboloid, as developed in [T51I16J . After that, we perform the triple sum estimates 
of the trilinear form using the individual decomposed estimates from Section [2j 

The generalized Littlewood-Paley projections are given by 

P,f{v) = / dv'2-^^{V{v - IL)) (v) f{v'), J > 0, 

JR" 

Q,/(z;)=>,/(t;)-P,_i/(iO, .?>1, 

where (p is a C°°, radial function supported on the unit ball of R"+^ chosen to 
satisfy various cancellation conditions (we refer to the discussion found in [TB] for 
exactly what is needed) and v = {v,^\v\'^) £ R"+^ for any w € R". Informally, 
Pj corresponds to the projection onto frequencies at most 2^ and Qj corresponds 
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to the projection onto frequencies comparable to 2-' (recall that the frequency 2^ 
corresponds to the scale 2~^ in physical space). We also define Qq = Pq- These are 
developed in J16j . In particular, we have Pjf{v) — >■ f{v) as j —)■ oo for all sufficiently 
smooth /. The principal reason for defining our Littlewood-Paley projections in 
this way is that the particular choice of paraboloid geometry allows us to control the 
associated square functions by our anisotropic norm. Specifically, in our previous 
papers [TB], it was established that 

(35) Y.'^''^'''^n dv\MQ,f{v)\ {vy^^^ <\f\l.,-,, (* = 0,1,2). 

This and more general results were proven in [121 Section 5 & (5.6)]. With that, 
we now establish Theorem [2] via the triple summation estimates in Section 13.11 



3.1. The main upper bound inequality. We will write / = YlT=ofj with the 
abbreviation Qjf = fj, and likewise for h. Now we expand the trilinear form 



oo oo 



oo oo 



(36) (Q(ff, /), /i) = E E (/j+'' Qah,) + E E (2;/^' '^^+') • 

1=1 j=0 1=0 j=0 

First consider the sum over I of the terms {fj+i, Qghj) for fixed j. We expand Q 
as a series by introducing the cutoff terms 2?^ and 2?^ from (fT6)) as follows: 

00 00 00 

E </^+'' 2^^^) = E E {^+(/.+'' '^.o - '^-(f,+i, h,)} 

l—l k— — oo l — l 

J 00 

(37) = E E{^+(/^+'''^^)-^-(/^+'>^)} 

(38) 



k— — oo l — l 

00 00 



+E E {'DXif,+i,h,)~viif,+i,h,)}. 

1=1 k=j + l 

Now the order of summation may be rearranged at will since the estimates we 
employ imply that the sum is absolutely convergent when h and / are Schwartz 
functions. Regarding the terms ([57]). the inequalities (PT|) and (P7)) dictate that 



One may now conclude by Cauchy- Schwartz that 

00 j 



j=Q fc= — 00 j=0 



<Cg2-'^ 



3=0 



|2 



E2"^i^.ii^, 



3=0 



<Cg2-''\f\N^,-.\h\Ns.-,. 



The comparison of the square function norm to the norm | • In".-/ is provided by 
?5|). This estimate may clearly be summed over / > 0. 
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To expand Q for the terms in (j36|l of the form {Q*gfj, hj+i) with the operators 
P^, T)^ and O* from ([19]) and (|20|) we use an analogous argument as follows: 



(39) =^a(/,>,+o+ E Y.^'DX-v':){f,,h,+i) 



1=0 k=-oo 1=0 



(40) +E E {'DX-v';){f,,h,+i). 

1=0 k=j + l 

The estimates ([M)) and (|27p are used to handle the terms in ([M)) involving the sum 
over fc < j just as for (|37p . except that that the roles of /i and / arc now reversed. 
To control the sum over O, (/j, ft.j+() in ([55]) . notice that 

oo oo c« oo 

/=o j=o 1=0 j=a 

which follows from Proposition [3] and the inequality 1 < 2'^"^ = 2-^'2'''(J+')2-'*-''. 
Again the Cauchy-Schwartz inequality and (f35|) yield the desired upper bound. 
Lastly, consider (p8| and (|40|). The terms (p8| are handled by (f32|) : we have 



Now there is decay of the norm as A: — )■ cso since 2s — z < 0, so that 

oo 

k=j+i 

Again Cauchy-Schwartz is applied to the sum over j. In this case 2'^^''~*)-' is written 
as 2('*~')J2'*(-'+'^2~'*'; the first factor goes with h, the second with /, and the third 
remains for the sum over I. Once again psp is employed. The desired bound for 
the trilinear term is completed by performing summation of the terms (|40p . The 
pattern of inequalities is the same, this time using p4p . In particular, one has 

1(2?'; - Pi-) (/.>.+0| < C,2(2-)^||Vr/,lL^^,J/i.+z|L^^,^. 
Again, as in Proposition [6l z ~ 1, 2 always satisfies 2s — z < 0, leading to 

oo 

Y^ |(p't-2?,')(/„/i,+0| <C,2(2-)^||V|7,lL^,,J/i,+dL^,,,- 

k=j+l 

The same Cauchy-Schwartz estimate is used for the sum over j; there is exponential 
decay allowing the sum over I to be estimated. The end result is Theorem [5J 

4. The coercive lower bound inequality 

The goal of this section is to prove PTj) in Theorem[TJ The main idea of the proof 
can be stated as follows: the principal analytical difference between the left- and 
right-hand sides of (fTTj) is that (for fixed ?;*) the variables v and v' are constrained 
relative to one another on the left-hand side, while on the right-hand side they 
essentially are not (they need only satisfy a distance inequality). To "regularize" the 
left-hand side and remove the constraint, we will exploit the fact that, for different 
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values of w*, the constraint (namely, the sphere) between v' and v changes. We will 
employ an elementary but completely novel convolution-type argument which will 
allow us to exploit the changing constraint. Loosely speaking, because the spheres 
arc not static, the convolution-type operation will "smear out" the support of the 
integral and give something analogous to the right-hand side. 

4.1. Regularization. As is customary, the proof proceeds by a careful dyadic 
decomposition. For any integer fc, consider the set flk — f2fc(w,w',w*) given by 

ilk = {iv,v',v^) I \v-v'\ < 2-^ and {2v' - v ~ v^,v - v^) > O}. 

Furthermore the condition (2u' — v — v^^^v — v<t) > is needed because of the sup- 
port condition in ([3|). We also use the quadratic functional given by 

h{f) = / dv, f dv f da,, (/' - /)2 In, g. \v - ^;,|("-l)+^+2^ 

JBn. Jm" "'Si;. 

By the integral da,', we mean /g„_i da (j){v') = L^, dayi(j){v'), where the left-hand 
side is exactly as in ([2]). In particular, note that E>^ is the sphere 

(41) §1^ = {w e M" \0={w-v,w-v,)}, 

which is the unique sphere for which v and w* are antipodal, having center ^"t^* 
and radius ^"~"*^ . With ([T5|) and ([7]) it is easy to check that 

oo 

(42) Ng{f)> Y, 2^«"-i)+2«)4.(/). 

k—~oo 

This follows because 

Kcos.)>(f^" 

\ |f — u' 

by virtue of ([3]) and ([2]); with this inequality, (j42|) follows from the observation 

CJO 

fc— — oo 

since the sum on the left-hand side is a geometric series which terminates at some 
maximal k satisfying 2^^™"^ sa |w — u'|~^ (and the sum is comparable to the value of 
the largest term). The goal will be to estimate the terms on the right side of (jH]) 
by something more directly comparable to our semi- norm ([6]). To that end, let 

z«fc(«,«*) = / da„,lo, <2-^("-i)|z;-«,r("-i), 

(the estimate from above follows because the angle to the pole is comparable to 
|w — w'||w — w, 1^^). Extending our convention by defining ^^ = g{v^,) and f = /(«'), 
Fubini's theorem guarantees that the quantity Ikif) Jg dv^^g,: is equal to both of 
the following integrals: 

I dvl duJ da^' I dv, I da,,{f - ff ^";^!!'°. g^-g, \v - ^;*|("-l)+^+2^ 

JK" J Br Js^^ J Br JsZ, Wk(V,V,) 

dvf dvj da^' I dv, I da,,(f' - ff ^";^"'°. g,-g, |« - TJ,|("-l)+^+2^ 
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where flk ~ ^kiv,v',v*)- We now bound Ik{f) Jg dv^g^^ below by an integral 
whose integrand is the average of the two integrands above. The elementary in- 
equality 

(43) if - If + a' - ff >\{y - n\ 

leads to the somewhat less elementary observation that 



I,, _ „ |("-l)+7+2s , I ^^ |(ri-l)+7+2s 

(44) (r-/)^ '" "*i _, — + (7'-/)2^Lj^iL — ^ — 

i\v-v,\\v~v,\)-i^-^) ^■' ■' ' ' 

where wc have also used the estimate of w^ from above. 

Since the quadratic difference now involves only v' and v' , the next step is to 
perform a change of variables simultaneously for both v' and v' (motived by the 
Carlcman-typc representations) so that, when wc integrate the right-hand side of 
(Pl)l over all the relevant variables, we may treat v' and v' as being unconstrained 
at the price of requiring v to simultaneously satisfy two constraints. Let H be any 
(Borel) measurable function oi v,v' and v' . Then for any fixed v* and ?7*, with the 
notation §^^ from (gT]) we have that 



(45) I dv I dajj' / da^, H{v,v',v') 



^ I dv' I d^' I d..— ^("'"''^') 



where 



"^ L/ „, 121777/ 777 |2 /„./ „. —I — \2 



(46) D = \v' - v,\^\v' -v,\^ - {v' - v,,v'^v 



*/ J 



and E^, '^* is the co-plane (by which we mean an affine subspace of codimension 2) 
of points V e M" satisfying the constraints 

(47) <,';:: = {t; e K" : = {v - v' ,v' ~ v,) , = {v ^v' ,v' - tJ,) }, 

and, as always, diTy is the Lebesgue measure on this co-plane. Now consider the 
quantity Kk — i^fc(u', u», w', w*) given by 

^''^ ""' = ^^4'::: — |.-..i-k-^*|-^ — '"^'^^- 

By the co-plane change of variables (|45|) , we may conclude that 

(49) hif) f dv,g, > f dv, g, f dv,g, f dv' f dv'Kk{f~l'f. 

Jbr Jbr Jbr Jr" Js." 

The goal at this point is, of course, to obtain a favorable estimate for Kk in terms 
of the geometry of the various variables of integration. 
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4.2. Estimation of (|48| . Regarding the function D from (j46|) . it is a useful fact 
to note tliat D is a Gram determinant, and its value is equal to 4 times the square 
of the area of the triangle with vertices 0, v' — w* and v' — w*. Consequently we 
have two alternate formulas for D which will also be useful: 

D = |w' — w*P|u' — u' + u, — w*P — (w' — w' + w, — u*,u' — u*) 

I—/ — \2\ I — /,— |2 // — 'i— — ' — \2 

~ \v — W,| \V —V +w*— W*| — (w — W +w*— W*,W — V^,) . 

For the moment, let us make a series of assumptions (to be examined later). To 
that end, we will consider w', u,, v', TJ* as fixed and satisfying the inequalities 

(50) |w*-w*| > (5, 

(51) D > smv'^v,\^ + \v'-v,\^), 

for some fixed constant 6 > 0. Now ((TT|) in Theorem [T] will be proved through a 
series of propositions estimating the various features of pS)) . 

Proposition 7. Assume that (|50p and (|5ip hold for some quadruple («',«*,«',«*) 
wii/i ^xed (5. T/ien /or a// fc > ko, with fco = ko{n, R, S) > 0, we have 

(52) /^,>2'=("-i)(.r+^^+^/ , d-. lo. In,- 

Proof. The inequality ([5^ rests on a variety of length inequalities which are a 
consequence of the assumptions ([50)) and ([5T|). First, assuming fc > 0, on Q^ n fifc 
we have that \v' — v'\ <2. By the lower bound for D, this implies that 

SWv' - v,\^ + \v' - w,n < \v' -v,\^\v' -v' + v, - v,\^ < {2R + 2)^\v'-v,\^, 

(where the bound from above for D is obtained by dropping the negative inner 
product squared in the third representation of D mentioned previously) so that 
\v' -v^\ < S~^{2R + 2)\v' -v^\. By symmetry, \v' -v^\ < 6-\2R + 2)\v' - v^\ as 
well. As long as k is large enough that 2^'^' < |, then, we also have that 

\v' — w,| + \v' — w,| > |w* — w,| — |w' — v'\ > -. 

In particular, we may conclude that |w'— w* | « \v'—v^, \ > 1 with constants depending 
only on S and R. So again, if k is large enough (depending on S and R) this means 
that any v in 51/; n fJ^ n E'^,'^* will satisfy \v — w, | w |u' — u, | « |ii — w, | « (v') 
(because in this case all of the lengths are bounded from below on (|T7|) . |u — u'| is 
bounded from above, and |w, |, |u*| are bounded above as well). Thus 

uniformly (where the constant depends on the dimension, i?, and 5 only). Finally, 
using the bound D < {2R + 2)^111' — w*p (obtained by keeping only the product of 
squared lengths and neglecting the squared inner product) gives (|52[) . D 



Lemma 8. Under the assumptions (j50p and (|5ip for some fixed quadruple of points 
{v' , u*, u',TJ*), there is a constant e > and a ko such that the additional constraint 
d{v' ,v') < 2~^£ for any k > ko implies that the set 

[v e El',f^ I \v - v'\ < 2-^ and \v - v'\ < 2-^= } , 
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is nonempty and contains a Euclidean hall of radius |2~ . 

Proof. The objective at this point is to show that, when v' and v' are close in 
the anisotropic sense, then the co-plane above passes near to these points in the 
isotropic sense. To that end, consider the temporary definitions 

dcf / 1 dcf / 

u = V — w*, and u = v — w*. 

Note that (pS)) is given simply by the expression D = |Mp|up — (w,u) . It is a 
straightforward exercise in linear algebra to verify that the vector w given by 

(53) w = v' + ^^' ~j''^^ (- {u,u) u + l^pn) , 

must lie in the co-plane E^, '^* from (|47|) (in fact, w is the point in the co-plane of 
minimal Euclidean distance to v'). The distance of this point w to v' is given by 



I '|2 

\W — V \ 



D 



{{u^uf \u\^ - 2 (u, w) \u\^ {u,u) + \u\^\uf\u\'- 



I -\2 
V ,U) 



D 

Expanding this inner product gives 

(w' - v ,v' - tJ*) = x(|wT ^ l^^'l^ ~ 1^' ~ "^T) ^ ^' - v ,v*) ; 
consequently, if d{v' ,v') < 1 then it follows that 

I- - -f < '"''^^*^y'^^^^' < {d{v',v')r, 

where the final constant depends on 5 and R. 

The crucial point is that if d{v',v') < £2^^ for e = e{n,6,R) > sufficiently 
small, then w will be within |2~'^ of v': \w — v'\ < |2~'"'. We may also choose e so 
that \v' — v'\ < |2^'°, whence the triangle inequality guarantees \v' — w\ < j2~'^. 
Another application of the triangle inequality shows that any v g E'^, '.^ * satisfying 
|u — ifl < |2~'^ will be forced to satisfy both \v — v'\ < 2^^ and |w — u'| < 2^'^. D 

Proposition 9. Under the same hypotheses as Lemma\^ we have the lower bound 

-f^fc ?C 2 {V ) i-d{v\v')<e2-i', 

uniformly for all k> k^, where ko ^ kf){n, S, R) > and e = e(n, S, R) > 0. 

Proof. When d{v',v') < e2~^, then Lemma [8] identifies a Euclidean ball of radius 
comparable to 2~^ which lies at a distance strictly less than 2~^ to both the points 
v' and v' . This is nearly sufficient to assert that this ball lies in the intersection 
rjfe n ilfc, but there is a second constraint to satisfy: namely, the positivity of the 
relevant inner products. For any v in this ball, on (|47p we have 

{2v' — V — v^.,v — v^,) ~ {v' — v^, — {v — v),v' — w* -|- (u — v)) 

I ' |2 I f\2 

= |w — w*| — \v — V \ . 

Again as long as k is sufficiently large, this quantity will be positive because jw'— w* p 
is bounded below. By symmetry, {2v' — w — tJ*, w — w*) will be positive in this case 
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fe 



as well. Consequently, for any fixed v',v' with d{v' ,v') < s2 , the ball of points 
V € E^, '^* identified in Lemma [5] will be contained in O^ n U,k '■ 

Here w is defined in ([5^ . Since the measure of this set in E^^/!^' is comparable to 
2-("-2)fe ^j|- jg codimension 2 in R"), we have arrived at the conclusion 

7^ > r,fe(n-l) /„/\7+2s+l r)-k(n-2)i 

uniformly for all k > ko, which is exactly the desired lower bound for (|3H]). D 

In light of (|49p and Proposition IH] we may now establish a temporary version 
of pT|) . In particular, by (|49p and Proposition IH] we have that, for all k > kg, 
Ikif)\g\L^{BR) is bounded uniformly below by 

(54) 2^= / du' / dv' U'-Tf («r^''^'ldK.^')<e2- A^,W,^'), 

with 

A'f = dv^ / dv^ g*g^lG- 



9 

^ Bu ^ Br, 



By ([42]) . one may estimate Ng{f) |5|li{Sh) from below by multiplying ([54]) by 
2(("-i)+2s)fe and summing over k. Interchanging summation and integration of the 
terms (j54p . we will next employ the inequality 



oo 



(55) J2 2'<"+''^ldK,^')<e2- > div',vT''-''Uv',-')<2-^os, 



k — kn 



which follows because, as before, the sum on the left is a finite geometric series (for 
any pair v' ,v' , the characteristic functions arc zero if 2~'"'e > d{v',v')). The sum 
of a (nontrivial) geometric series is comparable to its largest term, which, in this 
case, satisfies 2"'^™"^ « d{v',v'). Thus by ((55|) it follows that 



Tt's 



(56) N,{f) |.gk.(B«) > f dJj'f d.- j^. {|L i-y^"'^' l.K.^')<2-o. Af . 

JR" JR" "l^ 1 ^ J 

The proof of pTjl would be complete if we had a uniform positive lower bound for 
A'^{v',v'), and if we could replace d{v' ,v') < 2~^°e with d{v',v') < 1 in the above. 
Closing both of these issues is the content of the remainder of this section. 

4.3. Simplifying the geometry. The next task at hand is to gain a better geo- 
metric understanding of (|5l]) . Wc seek to replace this somewhat obscure condition 
with a more intuitive one. Specifically, we will replace (|50p and (|5ip with: 



(57) \v'-v4 > (5>0, 

(58) \v' -v,\^\v' -v,\^-{v' -v*,v' -v,f > S^\v'~v^\^. 

These new conditions have simple geometric interpretations: the distance from v' 
to w* is at least 5, and TJ* lies outside a tube of radius 6 around the line through v' 
and w*. 

Note also that ((58)) itself implies |w' — iJ*| > i5 as well (which may be seen by 
neglecting the squared inner product on the left-hand side). Likewise, the equality 

\v' — w,|^|w' — w*P — (w' — w*,w' — w») = |w* — tJ»P|u' — w*P — (w* — «*, w' — w») , 
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guarantees that |z;* — w*| > (5 as well. Now since \v' — TJ*| and \v' — v^\ are bounded 
below and the distances |i;* — TJ»| and \v' — v'\ are bounded above, it follows that 

jw' ^ TJh.| w |TJ' — TJ, I sa |ii' — w*| w |u' — u*| w {v') , 

when d{v' ,v') <^ 1; all of the implicit constants will depend on S and R. 

Proposition 10. The inequalities (|57p and (|58|) imply assumptions (|50p and (j5ip 
('wif/i a different S > 0) provided d[v' ,v') < e for some small e = e((5, i?). 

Proof. It remains only to verify ([5T|) . Note that ([55)1 coincides with ([51]) when 
v' = v' . It thus suffices to estimate the difference when going from v' to TJ' in the 
appropriate places. First, the change in the squared length \v' — v' + v^ — w, p is 



\v' -v' + v*-v4^ - !«*- w*n < \v' -v'\^ + 2\v' -v'\\v^ 



w* 



< 



V 



where we use the inequalities \v' — v'\ < 1 and |w* — w*] < 2R. Next, we observe 
the equality 

(f ' — U', w' — w») = (w' — if , v') — (v' — if , u*) 

^l{\v'\' -\v'\^ + \v' ~vf)- {v' ~v',v^) , 
meaning that | {v' — v',v' — v^) \ < d{v' ,v'). In particular, this implies 

I / ' — ' I — ' \ 2 /— / 

I ( W — W +IJ* — W*,W — W, ) — (i;, — W,,W — W,^ 

= |(w' — w',u' — w«)| Kw' — zJ' + zJ* — w*,w' — w*) + (u* — w*,i;' — «*)! 
Since lu' — i'*| > (5, it follows that the difference of the left-hand sides of (IST 



and (|58p are bounded in magnitude by a uniform constant times d{v' ,v')\v' — w*p 
(because |w' — w*| is bounded below, the squared length dominates the length itself). 
Assuming that d{v' ,v') <^ 1 and using |t;' — u*] « |tJ' — tJ,| establishes ([ST]). D 

Proposition 11. Suppose that ^ holds for some Cg > 0. Then 
Af= I dv, I dtJ, 1g 5, g, > Cl 

J Br ^ Br 

Proof. By Proposition [101 it suffices to estimate the integral A'? uniformly from 
below for all v' with the set G replaced by the set G" determined by the constraints 
([57]) and ([58]) . Clearly for each fixed w* with |t;* — w'| > 5 the integral over u* is 
exactly the integral of g over the ball Bji minus a tube of radius 5 determined by 
v' and Vst- Thus the Fubini theorem dictates that 

/ dv^ j dv^ 1q' g^g^ > Cg dv^, g^. 

J Br J Br JBr\Bs{v') 

To conclude, we observe that Bn minus the ball of radius S around v' is a strictly 
larger set than the ball minus any tube of radius S whose axis passes through v' . D 

At this point, with Proposition II ![ we may estimate A^ and consequently use 
([56[) to deduce pT|) aside from the limitation that the distance d{v' ,v') would be 
constrained to be less than or equal to e > rather than the larger bound 1 
appearing in the definition of N'^''^ . This turns out to be a minor difference which 
is fixed in the next subsection. 
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4.4. The completion of ([TT]) . For completeness, we eliminate the parameter 
2"'^°£ appearing in (|56|) and ([55]) . We will use the estimate (|43|) (with the roles of 
/ and / reversed) combined with the trick of adding extra variables of integration 
and bounding below. First we add extra variables of integration: 



\2 I /\7+2s+l 



Ld(u',t7')<e 



dv I dv' I du'if'-7'r{v'y^'-'+' 



J-d(v.v')<e 

-''^^''^'^^^ \B,iv')\ 



dv / dv' / dv'if -ff (v) 



7+2S+1 ^d{v',v')<e . 



<-d{v,v')<e- 



\Beiv')\ 

Here Be{v') is the anisotropic ball of radius e centered at v' and |i3j(TJ')| is its 
measure. Now we use (|43)) (together with the properties {v) « {v') and |i3e(w')| « 
|i3e(w')|, which hold uniformly when e < 1) to conclude that 

Jeif) + Jeif) > f dv f dv'if- f')^{v'y+''+'H{v,v'), 

where 



H{v,v') = -. / dv' ld{v.v')<e^d{v'.v')<e 



Now the equalities 





I' 4 


v' 










— v 






2 


v + v' 


1 


|7)P - 




2 




\v\ 





V + v' 



2' ' 2' ' 4' ' 

allow for the following estimate of the anisotropic distance: 

J,' v + v'Y \v-v'\^ 1 ( ( 1, „ 1, ,,2\ 1, ,,2 



1l,,2 , 1l/,2\ 1 



< 



d(zM^_ 1 

4 2 V 2'"' 2'^ ' y 4 

d(w,w')^ d(i;,w')^ d{v,v'f 

4 "^ 8 '^ 64 



i;-w'|2 + 



1 1 
416 



\v-v'\ 



/|4 



If d(w,w') < 1, then it follows that d ( «, '"~^" 1 < |(i(w,w'). By symmetry, it also 

follows that (i(w',^^^i^) < |d(f,v') as well. In particular, when d{v,v') < ^e, 
the intersection of d{v,v') < e and d{v' ,v') < e will contain the set of points v' for 
which d{v' , ^^^^j^) < ^s. Therefore, we have 



H{v,v')> 



\Be{v')\ 



d{v,v')<i§-e- 



Since e < 1, the ratio of the volumes above are uniformly bounded below (as both 
arc comparable to e" {v')~ ). We have therefore established that Jeif) > cJio^{f). 
for a constant c > which is uniform in e. Iterating this inequality, it follows that, 
for any fixed constant p > 

(59) Jeif) > CpJepif) 



26 P. T. GRESSMAN AND R. M. STRAIN 

provided only that pe < 1. Now by ([56]) and Proposition [TTJ we have 

oo 

But by ([59]). we have 

Je2-''if) > C2ko^-iJ2-k+ko, 

uniformly provided that fc > fcp, where the precise value of C2fcoe-i is irrelevant since 
it is ultimately determined only by the dimension, R, and 5. Thus we have that 

oo oo 

NM) \9\lhb,) ^ Cl Y. 2'^("+2^) J,2-.(/) > Cl Y 2''"'*-'^^ J2-+'=o(/). 

Shifting the index of summation on the right-hand side down by /cq and using the 
estimate ([55)1 with e = 1 clearly establishes that the right-hand side is uniformly 
bounded below by C^ |/|?v^ . This is exactly ([TT|) and therefore the proof of Theo- 
rem [T] is complete. 

5. Entropy production estimate 

In this section we prove Theorem [31 which will follow from Theorem [TJ We use 
the entropy production functional ([5]). From, for instance, [Il[ll], we split ([S]) as 

D{gJ)^SigJ)+TigJ), 

where 

S{g, /) "=' [ dv [ dv, [ da B{v -v,,o) g, (f log I + f - f) , 
T{gJ)'^^ - j dvt. g* I dv I da B{v-Vt,<7) (f ~ f) . 

JR" JR" JS"-i 

This follows easily from the prc-post coUisional change of variables. From the 
elementary inequality a log j — a + b > {\/a — Vb)"^, see e.g. [l], we obtain that 

S{gJ)> f dv f dv, f daB{v~v,,a)g, (^ - ^)\ 

which immediately implies S{g, f) > Ng{y/J) with Ng{^/f) from ([T]). Now Theorem 
[3] follows easily from Theorem [T] and the following estimate for Tg{f). Notice that 
with the cancellation lemma [T], using the arguments as in Section [131 with ([TJ and 
([Till , it is easy to see that |Tg(/)| < C"Cg|/|ii . This establishes Theorem[21 

Appendix: The dual formulation and other representations 

Finally, wc derive the dual formulation from PT]) and (1M| . For this we will use 
the Carleman-type representation which can be found for instance in |161 Appen- 
dix] . Notice also that [TH Appendix C] gives a proof of Carleman-type representa- 
tions. The functions 5 and $ below are given by Q and (HJ. However, to derive 
the dual formulation it suffices to suppose that both of these functions are smooth. 
The general expressions can then be deduced by approximation. Furthermore, we 
derive the "co-plane change of variables" from ps]) . 
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Dual Representation. We initially suppose that /g„_i dcr |fo((fc, cr))| < oo and 
that the kernel b has mean zero, i.e., /g„_i da b{{k,a)) = 0. Then after the pre- 
post change of variables (w,w*) — > («',«!) we can express (P) as 

{f.Qgh)= [ dv [ dv, [ damv~v,\)b{{k,a))gJih'-h) 



dv / du, / da ^{\v - v^\)b{{k,(T))g^fh'. 

Jr" JS"-i 

Now with the Carleman representation, e.g. |16[ Appendix], we have 

hi/ ^~^* 2v' — v—v^ \ \ 

(/, Q,h) = 2-1 / dv, f dv' f dn, mv - v4) )^J^^'!'^'~^'7:li^ 9*.fh'. 



E^l \v' - V*| \v - W*|" 2 

The definitions of these notations, E^ and dn^, were given previously in the para- 
graph containing ([T7)) . Furthermore from the identity (P5|) we observe that 

, , , v — Vf 2v' — V — v^, \\ Iw' — w, 1"^-^ 
diTy b ' ' ^ ' 



\v — w, I ' \2v' — V — v,\ / J jw — i;*P"-2 



da I r"-2 dr b f ''. ''* ' '^^ 1 . '", . "*' ,.. , = 0, 



v|w'-W,p+r2y (^2 _^ 1^/ _^^|2)n-l 



by a change of variables to polar coordinates since J_^dt b{t){l — t^) 2 =0 
(following from the cancellation condition on S"^^) and 



d 
dr 



\v' — w*P — r^ 



\v' — u* r + r^ 



-4r\v' - w*p 
(r2 + |v'-w,|2)2' 



71-3 



'-- |2_.2\^\ - (2r|t;'-'y,|)"-3 



/ (r2 + |i;' - i;,|2)"-3- 

In particular, this implies with B from (jlSp that 

a7r„ B -— ^ f-^ — g, f h = 0. 

We subtract this expression from the Carleman representation just written for 
(/, Qgh), and use the kernel defined in ([T5|) to see that (|17p holds with p^ . 

The claim is now that this representation (|17p holds even when the mean value 
of the singular kernel b{{k, a)) from ([3]) is not zero. To see this claim, suppose that 
b is integrablc but without mean zero. Then define 

b,{t)''^'b{t)-i[,_,^^{t)l ^dtb{t){i-f)^ ^' dt{i~e)^^ . 

As a function on E>"~^, b^ will clearly have a vanishing integral. However, given 
arbitrary /, g and h which are Schwartz functions, it is not hard to sec that 

Above Qi is the operator Qg formed with b^{t) in place of b{t). This convergence 
holds because cancellation guarantees that the integrand vanishes on the set defined 
by (fcjCT) = 1. Moreover, an additional cutoff argument shows that the equality 
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also holds provided that b{t) satisfies ([3]); the higher-order cancellation is preserved 
because i" _J* | possesses radial symmetry in v — v' . 

The "dual representation" ((T7)) deserves its name because if one defines 

Qgh{v) ^ / du, / d<7 B 5, {h' — h) , 



S"-i 



dof 



$(u' - v^)\v' 



then 

(60) {Q{9J),h) = {f,Q,h) = {Qlf,h). 

Note that the last inner product above represents an integration over dv' whereas 
the first two inner products above represent integrations over dv. 

The co-plane identity. The goal of this section is to prove the co-plane change 
of variables from (|45p . It suffices to prove the formula (|45|) for continuous, com- 
pactly supported H. The proof follows rather directly from the coarea formula for 
codimension 2, see e.g. |13] . For fixed v, for example, the coarea formula gives 

day>H= lim TT— ] - dv' H\v - v^.\luy,_yy_yAi<±, 

since the sphere S^^, as a function of v' , is precisely the zero set of the Lipschitz 
function {v' — v,v' — v^) as in (|4ip . The quantity |z; — t;^,| inside the integrand is 
the magnitude of the gradient of {v' — v,v' ~ v<,) on its zero set (which is required 
by the coarea formula) and the external term 2"~^|w — u*|~'"~^' is a normalization 
factor. Thus the left side of (|l5|) may be realized as a limit as e — )• 0''' of 

Now use the Fubini theorem to evaluate the integral with respect to v first. In this 
case the limit is supported on the intersection of the two zero sets, which is exactly 
the co-plane E^, '^* . The 2-dimensional Jacobian of these two constraint functions 
(now with respect to w) is exactly D^ from (|15)) . so that we have 



hm ^r / dv- 



0+ e^ ./»" \v — w,|" 2|w 



diiv- 



'El/l* j-y — w*|" 2|w — i'»|" 2-D2 
Substituting this back into the integral over v' and v' establishes (|45|) . 
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